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We show that quantum feedback control can be used as a quantum error correction process for 
errors induced by weak continuous measurement. In particular, when the error model is restricted 
to one, perfectly measured, error channel per physical qubit, quantum feedback can act to perfectly 
protect a stabilizer codespace. Using the stabilizer formalism we derive an explicit scheme, involving 
feedback and an additional constant Hamiltonian, to protect an (n — l)-qubit logical state encoded 
in n physical qubits. This works for both Poisson (jump) and white-noise (diffusion) measurement 
processes. In addition, universal quantum computation is possible in this scheme. As an example, 
we show that detected-spontaneous emission error correction with a driving Hamiltonian can greatly 
reduce the amount of redundancy required to protect a state from that which has been previously 
postulated [e.g., Alber et al, Phys. Rev. Lett. 86, 4402 (2001)]. 
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I. INTRODUCTION 

Many of the applications of quantum information sci- 
ence, such as quantum computation P, and quan- 
tum cryptography rely on preserving the coherence 
of quantum states. However, these states are typically 
short-lived because of unavoidable interactions with the 
environment. Combatting this decoherence has been the 
subject of much study. 

Two important tools that have been developed for this 
task are quantum error correction 0, |E IE and quan- 
tum feedback |^ |^, ^| . In the usual protocol for quan- 
tum error correction, projective measurements are per- 
formed to acquire an error syndrome. A unitary opera- 
tion chosen based on the results of the projective mea- 
surements is then applied to correct for the error. Quan- 
tum feedback control, on the other hand, uses the tools 
of continuous measurements and Hamiltonian feedback. 
The parameter to be controlled is typically the strength 
of the feedback Hamiltonian, which is conditioned on the 
result of the continuous measurements. 

Quantum error correction and quantum feedback both 
rely on performing operations that are conditioned on 
the result of some measurement on the system, which 
suggests that exploring the links between these two tech- 
niques adds to our understanding of both processes, and 
may lead to insights into future protocols and experimen- 
tal implementations. In particular, this work provides an 
alternate avenue for examining the situation considered 
in ^1,12, 13,2^] of correcting for a specific error process, 
such as spontaneous emission, at the expense of correct- 
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ing fewer general errors. Practically, as these authors 
point out, it makes sense to pursue the tradeoff between 
general correction ability and redudancy of coding, as 
smaller codes are more likely to be in the range of what 
can be experimentally realized in the near future. We 
shall see that combining the pictures of quantum feed- 
back and error correction provides a convenient frame- 
work in which to investigate this situation. 

An additional motivation for considering the union of 
these techniques, as in |l5j |. is to examine what is possi- 
ble with different physical tools: in particular, continuous 
measurements and Hamiltonians instead of the projective 
measurements and fast unitary gates generally assumed 
by discrete quantum error correction. Continuous error 
correction might well be useful even in a scenario in which 
near-projective measurements are possible (e.g., ion traps 
^3 and superconducting qubits [l3|); it could be modi- 
fied to provide bounds on how strong interactions in such 
systems would have to be to perform operations such as 
error correction and stay within a certain error threshold. 

Ref. 13 presupposed that classical processing of cur- 
rents could be done arbitrarily quickly, so the feedback 
was allowed to be an extremely complicated function of 
the entire measurement record. This can be modeled only 
by numerical simulations. In this paper, by contrast, we 
will restrict our feedback to be directly proportional to 
measured currents, thus removing any need for classi- 
cal post-processing. In the Markovian limit, this allows 
an analytical treatment. This simplification is possible 
because in this paper we assume that the errors are de- 
tected. That is, the experimenter knows precisely what 
sort of error has occurred because the environment that 
caused the errors is being continuously measured. Since 
the environment is thus acting as part of the measure- 
ment apparatus, the errors it produces could be consid- 
ered measurement-induced errors. 

There are a number of implementations in which 
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measurement-induced errors of this sort may be signif- 
icant. In the efficient hnear optics scheme of Knill et al. 
[Tsj l. gates are implemented by nondeterministic telepor- 
tation. Failure of the teleportation corresponds to a gate 
error in which one of the qubits is measured in the com- 
putational basis with known result. In a number of solid 
state schemes, the readout device is always present and 
might make an accidental measurement of a qubit, even if 
the readout apparatus is in a quiescent state. An example 
is the use of RF single electron transistors to readout a 
charge transfer event in the Kane proposal. Such a mea- 
surement is modelled as a weak continuous measurement 
[Tgl . While one supposes that the SET is biased in its low 
conductance state during qubit processing, it is useful to 
know that even if the device does accidentally make a 
measurement, the resulting error can be corrected. 

In this paper, we show that for certain error models 
and codes, Markovian feedback plus an additional con- 
stant Hamiltonian (a "driving Hamiltonian" ) can pro- 
tect an unknown quantum state encoded in a particular 
codespace. Using the stabilizer formalism, we show that 
if there is one sort of error per physical qubit, and the 
error is detected perfectly, then it is always possible to 
store n — 1 logical qubits in n physical qubits. This works 
whether the detector record consists of discrete spikes 
(Poisson noise) or a continuous current (white noise). 
This suggest that if the dominant decoherence process 
can be monitored, then using that information to control 
the system Hamiltonian may be the key to preventing 
such decoherence (see also the example in [20|'). 

As a salient application of this formalism, we con- 
sider the special case of spontaneous emission. Sta- 
bilizing states against spontaneous emission by using 
error-correcting codes has been studied by several groups 
[iTl IT^ IT^ Here we demonstrate that a simple 

n-qubit error-correcting code, Markovian quantum feed- 
back, and a driving Hamiltonian, is sufhcient to correct 
spontaneous emissions on rt — 1 qubits. The result of 
encoding n—1 logical qubits in n physical qubits has 
been recently independently derived by |2l| for the spe- 
cial case of spontaneous emission; however, our scheme 
differs in a number of respects. We also show that spon- 
taneous emission error correction by feedback can be in- 
corporated within the framework of canonical quantum 
error correction, which can correct arbitrary errors. 

The paper is organized as follows. We review some 
useful results in quantum error correction and quantum 
feedback theory in SecE] In Sec. Illll we present the ex- 
ample of detected spontaneous-emission errors, first for 2 
qubits and then for n qubits. In Sec. II VI we generalize this 
for protecting an unknown state subject to any single- 
qubit measurements. We show how to perform universal 
quantum computation using our protocol in Sec. Sec. 
I VII concludes. 



II. BACKGROUND 
A. Quantum error correction: stabilizer formalism 

Quantum error correction has specifically been de- 
signed for protecting unknown quantum states 6, 7,. .2^ . 
An important class of quantum error-correcting codes are 
the stabilizer codes. An elegant and simple formalism 
exists for understanding these codes; in this paper we will 
restrict ourselves to this class of codes in order to take 
advantage of this formalism. 

In the remainder of this paper we will use the notation 
of "^l in which X, Y, and Z denote the Pauli matrices 
ax, cTy and az respectively, and juxtaposition denotes a 
tensor product; hence any element of the Pauli group 

P,,^{±l,±i}(g>{I,X,Y,Z}^" (2.1) 

may be denoted as a concatenation of letters (e.g., ZZI = 

A stabilizer code may be defined simply as follows: 
Consider a 2"-dimensional (n-qubit) Hilbert space and 
a subgroup of 2""*^ commuting Pauli operators S G P„. 
This group of operators is the stabilizer of the code; the 
codespace C{S) is the simultaneous -1-1 eigenspace of all 
the operators in S. It can be shown that if — / is not an 
element of 5, the subspace stabilized is non-trivial, and 
the dimension of C{S) is 2*^; hence, we regard this system 
as encoding k qubits in n. The generators of such a group 
are a subset of this group such that any element of the 
stabilizer can be described as a product of generators. 
It is not hard to show that n — k generators suffice to 
describe the stabilizer group S. 

When considering universal quantum computation it 
is also useful to define the normalizer of a code. Given a 
stabilizer group 5, the normalizer N{S) is the group of 
elements in P„ that commute with all the the elements 
of S, and it can be shown that the number of elements 
in N{S) is 2"+*=. 

Now, n + k generators suffice to describe N{S). Of 
these, n — k can be chosen to be the generators of S. It 
can be shown that the remaining 2k generators can be 
chosen to be the encoded operators Z^, X^, fj, = 1,2, k, 
where Z^^, X^^ denote the Pauli operators X and Z acting 
on encoded qubit ^, tensored with the identity acting on 
all other encoded qubits. These encoded operators act, 
as their name implies, to take states in C{S) to other 
states in C{S). 

The usual protocol for stabilizer codes, which will be 
modified in what follows, starts with measuring the stabi- 
lizer generators. This projection discretizes whatever er- 
ror has occurred into one of 2""*^ error syndromes labeled 
by the 2""'^ possible outcomes of the stabilizer generator 
measurements. The information given by the stabilizer 
measurements about what error syndrome has occurred 
is then used to apply a unitary recovery operator that 
returns the state to the codespace. 

In this paper we will use a modified version of this 
protocol. In particular, we will not measure stabilizer 
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elements. Instead, we will assume that a limited class of 
errors occurs on the system and that these errors are de- 
tectable: we know when an error has happened and what 
the error is. The correction back to the codespace can 
still be performed by a unitary recovery operator based 
on the information from the error measurement. Fig. ^ 
shows the difference between the conventional protocol 
and our modified protocol. 

Conventional protocol 
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FIG. 1: The top diagram shows the conventional stabilizer 
error correction protocol. After the state is encoded, an error 
occurs through coupling with the environment. To correct 
this error, the encoded state is entangled with a meter in or- 
der to measure the stabilizer generators, and then feedback 
is applied on the basis of those measurements. The bottom 
diagram shows our modified protocol, in which the error and 
measurement steps are the same. To correct the error in this 
protocol, the environment qubits are measured, and we feed- 
back on the results of the environment measurement. 

In this paper, we will also consider operators of the 
form 



considering the dynamics of the measured system condi- 
tioned on the continuous measurement record; this pro- 
cess is referred to as unraveling. The reduced dynamics of 
a system subject to weak continuous measurement is de- 
scribed by a Markov master equation, which determines 
the dynamics of the system averaged over all possible 
measurement records. However, if the time-continuous 
measurement record (a classical stochastic process) is 
known, then it is possible to describe the conditional 
state of the measured system by a stochastic conditional 
evolution equation. A given master equation does not 
uniquely determine the conditional evolution equation, 
as there are many ways in which information about the 
system may be collected from the environment to which 
it is coupled as a result of the measurement. That is to 
say, a given master equation admits many unravelings. 

In this section we will introduce some of the results of 
this formalism; for more details see We will assume 
that the change in the state of the system over a time 
interval dt due to its interaction with the environment 
can be described by a single jump operator c. By this we 
mean that jumps are represented by a Kraus operator 
r^i = c^/dt, so that they occur with probability {c'cjdt. 
Normalization requires another Kraus operator, fio = 
1 — cJcdt/2 — iHdt, where H is Hermitian. Then the 
unconditional master equation without feedback is just 
the familiar Lindblad form 

dp = ilopflo + flipili — p 

= —i[H, p]dt + cpc^dt — —{c)cp + pc^c)dt 



-i[H, p]dt + V[c]pdt. 



(2.3) 



A bosonic example is given in |25j, while a fermionic ex- 
ample is given in p^ . 



T = Ti 



> Tr,. 



(2.2) 



where Ti is an arbitrary traceless one-qubit operator nor- 
malized such that its eigenvalues are {—1, 1}. Operators 
of this form are not generally Pauli-group stabilizers as 
presented in 0, as T is not in general a member of P„. 
However, because of the special form of T, T is equiva- 
lent to a Pauli operator up to conjugation by a unitary 
that is a product of one-qubit unitaries, i.e., there exists 
some U ~ Ui such that UTW is a member of P„. 

Therefore, choosing T as the sole stabilizer generator for 
a code is equivalent, up to conjugation by a unitary, to 
choosing a member of the Pauli group as the stabilizer 
generator. (Note that additional constraints are neces- 
sary if T is not the only stabilizer generator.) 



B. Quantum feedback 

Continuous quantum feedback can be defined, for the 
present purposes, as the process of monitoring a quantum 
system and using the continuous (in time) measurement 
record to control its dynamics. It can be analysed by 



1. Jump unravelings 

One way to unravel this master equation is to assume 
that the environment is measured so that the time of 
each jump event is determined. If the measured number 
of jumps up to time t is denoted N(t), then the increment 
dN{t) is defined by 



dN,{t)^ 
E[dN,{t)] 



dN,{t) 
{c^c)cdt. 



(2.4) 
(2.5) 



Here E[ ] defines a classical ensemble average, and the 
subscript c on the quantum average reminds us that the 
rate of the process at time t depends on the conditional 
state of the quantum system up to that time. That is 
to say, it depends on the state of the quantum system 
conditioned on the entire previous history of the current 
dN/dt. This conditional state is determined by a stochas- 
tic Schrodinger equation 



d\Mt)) 



dN,{t) 



dt 
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\Mt))(2.6) 



We will refer to this as a jump unraveling. If we 
average over the measurement record to form p{t) = 
E[\^c{t)) {i^c{t)\], it is easy to show using Eqs. (|2.4|) and 
(|2.5|l that p{t) obeys the unconditional master equation 
given in Ea. (|2.3|l . 

Now consider Markovian Hamiltonian feedback, linear 
in the current: 



(2.7) 



with V an Hermitian operator. Taking into account that 
the feedback must act after the measurement, it can be 
shown |8| that the feedback modifies the conditional evo- 
lution by changing the c in the numerator of the first 
term into e~*^c. Since likewise changing all of the other 
occurrences of c has no effect, the ensemble average be- 
haviour is the same as before, with c changed to e^'^c. 
That is to say, the feedback-modified master equation is 



p = -i[H,p\+V[e-''^ c\p. 



2. Diffusive unravelings 



(2.8) 



A very different unraveling may be defined by first not- 
ing that given some complex number 7 = |7|e*'^, we may 
make the transformation 



c 
H 



c -I- 7 

rr *l7L - 



(2.9) 



and obtain the same master equation. In the limit as 
I7I becomes very large, the rate of the Poisson process is 
dominated by the term |7p. In this case it may become 
impossible to monitor every jump process, and a better 
strategy is to approximate the Poisson stochastic process 
by a Gaussian white-noise process. 

For large 7, we can consider the system for a time St 
in which the system changes negligibly but the number 
of detections 6N{t) « Ijl^St is very large; then we can 
approximate SN{t) as [27] 



SN{t) 



|7|(e-**c + cte*^) 



St+\j\SWit), (2.10) 



where SW(t) is normally distributed with mean zero and 
variance St. 

We now define the stochastic measurement record as 
the current 



dQjt) 
dt 



lini 

7 — >oo 



SN{t) - \j\^St 
\l\5i 



(2.11) 



{e-"^c + e^'^c^l + dW{t)/dt. (2.12) 



Given this stochastic measurement record, we can deter- 
mine the conditional state of the quantum system by a 



stochastic Schrodinger equation analogous to Eq. (|2.6|l . 
The equivalence (in the ensemble average) to the master 
equation 1)2. 3|l is, in this case, easier to see by consid- 
ering pc — \'ipc){'4'c\, which obeys the stochastic master 
equation 

dpc{t) = -i[H, pc{t)]dt + V[e-"^c]pc{t)dt 

+ n[e-''^c]pc{t)dW{t). (2.13) 

In the above equations, the expectation (a)^ denotes 
tr (pcCi) , dW is a normally distributed infinitesimal ran- 
dom variable with mean zero and variance dt (a Wiener 
increment 2.8]), and 7i is a superoperator that takes a 
jump operator as an argument and acts on density ma- 
trices as 



Ti.[c]p = cp + pc^ — p tr [cp + pc'] 



(2.14) 



We thus have a different unraveling of the original mas- 
ter equation Ea. (|2.3|l . Because of the white noise in the 
stochastic master equation (|2.13ll we call this a diffusive 
unraveling. It applies, for example, when one performs a 
continuous weak homodyne measurement of a field c by 
first mixing it with a classical local oscillator in a beam- 
splitter and then measuring the output beams with pho- 
todetectors 0]. In that case the measurement process 
dQ{t) determines the observed photocurcnt. Another 
measurement model in which it may be appropriate to 
approximate a Poisson measurement process by a white- 
noise measurement process is the electronic point con- 
tact model for monitoring a single quantum dot[29ll30l|. 
In that case the form of the master equation itself de- 
termines a large background jump rate, rather than an 
imposed classical field prior to detection. 

Wc now consider Markovian feedback of the white- 
noise measurement record via a Hamiltonian, where the 
strength of the feedback is a linear function of the mea- 
surement current: 



Hfb{t) 



dQjt) 
dt 



(2.15) 



where is a Hermitian operator. It can be shown that 
the addition of such feedback leads to the conditioned 
master equation [8Ll3l| 



p = -i[{e">'c''F + e-"l'Fc)/2 + H,p] 



+V[e 



iF]p 



+dW{t)l-L[e~"t'c-iF]p. 



(2.16) 



In order to derive analytic results given such feedback, it 
is convenient to consider the average over many such evo- 
lution trajectories. Since the expectation value of dW is 
zero, averaging yields an unconditioned master equation 



-i[{e"^c''F + e^"^Fc) /2 + iJ, p] 



+V[e 



-i(f) 



iF]p 



(2.17) 



Note that these equations are only valid for perfect (unit- 
efficiency) detection; the correspondences between error 
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correction and feedback are more readily seen in this case, 
and we discuss the case of imperfect detection in Sec. 

MB 

These feedback equations are easily generalized in the 
following way: Given n qubits, denote a set of mea- 
surement operators by {ci, C2, • • • , c„}, where Cj acts 
on the jth qubit, and a set of feedback operators by 
{i^i, • • • , F„}, where the action of Fj is conditioned on 
the measurement of the jth qubit. Then the uncondi- 
tional master equation 12.17|l . for example, generalizes 
to 

n 

+V[e"''^Cj -iFj]p}. (2.18) 

III. EXAMPLE: SPONTANEOUS-EMISSION 
CORRECTION 

A particular example of a Poisson process error is spon- 
taneous emission, in which the jump operator is propor- 
tional to |0)(1|, so that the state simply decays from |1) 
to |0) at random times. Indeed, if the decay is observed 
(say by emitting a photon which is then detected), this 
may be regarded as a destructive measurement of the 
operator |1)(1|. 

Stabilizing states against the important decay process 
of spontaneous emission through application of error- 
correcting codes has been studied by several groups 
[ll, 12, 13, l||. In [13 Plenio, Vedral and Knight con- 
sidered the structure of quantum error correction codes 
and addressed the problem that spontaneous emission 
implies continuous evolution of the state even when no 
emission has occurred. They developed an eight-qubit 
code that both corrects one general error and corrects 
the no-emission evolution to arbitrary order. 

More recently, in several papers Alber et al. jT^ 
have addressed a somewhat more specific problem relat- 
ing to spontaneous emission from statistically indepen- 
dent resevoirs. In this formulation, the only errors pos- 
sible are spontaneous emission errors, and the time and 
position of a particular spontaneous emission is known. 
They showed that given these constraints, a reduction of 
the redundancy in ^3 was possible, and constructed a 
four-qubit code which corrects for one spontaneous emis- 
sion error. 

Here we show that for the case considered in , 
a very simple error correcting code consisting of just two 
qubits with feedback is sufficient to correct spontaneous 
emissions for a sing le logical qubit. A crucial difference 
from Refs. [isL IT^ is that we call for a constant driv- 
ing Hamiltonian in addition to the feedback Hamiltonian. 
Moreover, a simple code of n qubits, with the appropri- 
ate feedback and driving Hamiltonians, can encode n — 1 
qubits and correct for spontaneous emissions when the 
position (i.e. which qubit) and time of the jump are 
known. We also show that an equally effective protocol 



can be found for a diffusive unraveling of the spontaneous 
emission (as in homodyne detection). 

A. Two-qubit code: Jump unraveling 

The simplest system for which we can protect against 
detected spontaneous emissions is a system of two qubits. 
We consider the model in which the only decoherence 
process is due to spontaneous emission from statistically 
independent reservoirs. We will show that a simple code, 
used in conjunction with a driving Hamiltonian, protects 
the codespace when the time and location of a sponta- 
neous emission is known and a correcting unitary is ap- 
plied instantaneously; the codespace suffers no decoher- 
ence. 

The codewords of the code are given by the following: 

|0) ^ (|00) + |11»/V2 

|T) ^ (|01> + |10»/V2. (3.1) 

In the stabilizer notation, this is a stabilizer code with 
stabilizer generator XX. Both codewords are -t-1 eigen- 
states of XX. 

Following the presentation in Sec. Ill Bl the jump oper- 
ators for spontaneous emission of the jth qubit are 

Vlj — Kjdt{Xj — lYj) = -y/ Kjdtaj, (3.2) 

where 4k j is the decay rate for that qubit. In the absence 
of any feedback, the master equation is 

p = K,V[X, - iY,]p - i[H, p\. (3.3) 

If the emission is detected, such that the qubit j from 
which it originated is known, it is possible to correct back 
to the codespace without knowing the state. This is be- 
cause the code and error fulfill the necessary and suffi- 
cient conditions for appropriate recovery operations : 

(V-^l^t^lT/-,) = AbV. (3.4) 

Here E is the operator for the measurement (error) that 
has occurred and A^ is a constant. The states IV'm)' IV'i^) 
are the encoded states in Eq. (|3.1|l with — Sf_i,y. 

These conditions differ from the usual condition only by 
taking into account that we know a particular error E = 
rij has occurred. 

More explicitly, if a spontaneous emission on the first 
qubit occurs, |0) — *■ |01) and |T) — *■ |00), and similarly for 
spontaneous emission on the second qubit. Since these 
are orthogonal states, this fulfills the condition given in 
(|3.4() . so a unitary exists that will correct this sponta- 
neous emission error. One choice for the correcting uni- 
tary is 

Ui = {XI-ZX)/V2 

U2 = [IX - XZ)/V2. (3.5) 
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As pointed out in a further complication is the 
nontrivial evolution of the state in the time between spon- 
taneous emissions. From Sec. Ill Bl this is described by 
the measurement operator 

ilo = // (1 - (ki + K,2)dt) - KldtZI 

-K2dtIZ - iHdt. (3.6) 

The non-unitary part of this evolution can be corrected 
by assuming a driving Hamiltonian of the form 

H = -{kiYX + K2XY). (3.7) 

This result can easily be seen by plugging H3.7(l into (|3.6|l 
with a suitable rearrangement of terms: 

no = II{l-{Kl+K2)dt)-KldtZI{II-XX) 

- K2dtIZ{II - XX), (3.8) 

and since // — XX acts to annihilate the codespace, flo 
acts trivially on the codespace. 

We then have the following master equation for the 
evolution of the system: 

dp — UopUq — p + dt KjUjajpa!jUj , (3.9) 

i={l:2} 

where Uj is the recovery operator for a spontaneous emis- 
sion from qubit j. From Sec. Ill Bl these unitaries can be 
achieved by the feedback Hamiltonian 

where Nj{t) is the spontaneous emission count for qubit 
j, and Uj = eKp{—iVj) . Here, we can see from the simple 
form of 13.5|l that Vj can be chosen as proportional to Uj . 
Since UjajpOjUj acts as the identity on the codespace by 

definition, and since we have shown that flopilQ preserves 
the codespace, (|3.9f) must preserve the codespace. 

Such a code is optimal in the sense that it uses the 
smallest possible number of qubits required to perform 
the task of correcting a spontaneous emission error, as 
we know that the information stored in one unencoded 
qubit is destroyed by spontaneous emission. 



B. Two-qubit code: Diffusive unraveling 

A similar situation applies for feedback of a continuous 
measurement record with white noise, as from homodyne 
detection of the emission. We use the same codewords, 
and choose = — 7r/2 for the measurement. Then (|3.5() 
suggests using the following feedback operators: 

Fi = ^(XI-ZX) 

F2 = y^ilX^XZ). (3.11) 



If we use these feedback Hamiltonians with the same driv- 
ing Hamiltonian 1)3. 7|l as in the jump case, the resulting 
master equation is, using l|2.18|l . 

p = KiV[YI - iZX]p+ K2V[IY - iXZ]p (3.12) 

We can see that this master equation preserves the 
codespace, by again noting that YI — iZX — YI{II — 
XX), and similarly for lY-iXZ. The operator II-XX 
of course acts to annihilate the codespace. This insight 
will be used in the next section to derive a feedback pro- 
cedure for a more general measurement operator. 

C. Generalizations to n qubits 

We will now demonstrate a simple n-qubit code that 
corrects for spontaneous emission errors only, while en- 
coding n~l qubits. Both of the above calculations (jump 
and diffusion) generalize. The master equation is the 
same as (|3.3f) . but now the sum runs from 1 to n. Again 
we need only a single stabilizer generator, namely X®". 
The number of codewords is thus 2"~^, enabling n — 1 
logical qubits to be encoded. Since it uses only one phys- 
ical qubit in excess of the number of logical qubits, this 
is again obviously an optimal code. 

First, we consider the jump case. As in Sec. IIII Al 
a spontaneous emission jump fulfills the error-correction 
condition H3.4|l (see Sec. IIV Al below) . Therefore, there 
exists a unitary that will correct for the spontaneous- 
emission jump. Additionally, it is easy to see by analogy 
with that 

H = KjJ2 X^^-^YX"^"-^ (3.13) 
j 

protects against the nontrivial no-emission evolution. 
Therefore the codespace is protected. 

Next, for a diffusive unraveling, we again choose (f>j — 
— 7r/2, as in Sec. 1111 Bl The same driving Hamiltonian 
(|3.13|) is again required, and the feedback operators gen- 
eralize to 

= + x^^-^ZA:®"-^') . (3.14) 

The master equation becomes 

p = ^KjI?[/®J-ir/®"-^ (/®" - X®")]. (3.15) 

j 

These schemes with a driving Hamiltonian do not have 
the admittedly desirable property of the codes given in 
[1^ [13I IT^ that if there is a time delay between the 
occurrence of the error and the application of the cor- 
rection, the effective no-emission evolution does not lead 
to additional errors. Nevertheless, as pointed out in [T^ . 
the time delay for those codes must still be short so as 
to prevent two successive spontaneous emissions between 
correction; they numerically show that the fidelity decays 
roughly exponentially as a function of delay time. There- 
fore, we believe that this drawback of our protocol is not 
significant. 
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IV. ONE-QUBIT GENERAL MEASUREMENT 
OPERATORS 

The form of the above example strongly indicates that 
there is a nice generalization to be obtained by consider- 
ing stabilizer generators in more detail. In this section, 
we consider an arbitrary measurement operator operat- 
ing on each qubit. We find the condition that the sta- 
bilizers of the codespace must satisfy. We show that it 
is always possible to find an optimal codespace (that is, 
one with a single stabilizer group generator). We work 
out the case of diffusive feedback in detail and derive it 
as the limit of a jump process. 



A. General unraveling 

Different unravelings of the master equation 1)2. 3|l may 
be usefully parameterized by 7. In Sec. IIIBI we have 
seen that a simple jump unraveling has 7 = 0, while the 
diffusive unraveling is characterized by I7I — > 00. We will 
now address the question of when a unitary correction 
operator exists for arbitrary 7, i.e., when a measurement 
scheme with a given 7 works to correct the error. 

Consider a Hilbert space of n qubits with a stabilizer 
group {Si\. Let us consider a single jump operator c 
acting on a single qubit. We may then write c in terms 
of Hermitian operators A and B as 

e-"t'c = xI + A + iB (4.1) 
= xl + S-^ + ib-a (4.2) 

where % is a complex number, a and b are real vectors, 
and (7 = {X, Y, Z)^. 

We now use the standard condition (|3.4|) . where here 
we take E — c + j. Henceforth, 7 is to be understood 
as real and positive, since the relevant phase (j) has been 
taken into account in the definition H4.1|l . The relevant 
term is 

+Re(x + 7)A + Im(x + j)*iB + {a x b) ■ a) 
= i\x + l\^ + S^ + b^)I + D, (4.3) 

where D is Hermitian. 

Now we can use the familiar sufficient condition for 
a stabilizer code 0: the stabilizer should anticommute 
with the traceless part of E^E. This condition becomes 
explicitly 

Q^{S,D}. (4.4) 

As long as this is satisfied, there is some feedback unitary 
e~'^ which will correct the error. 

Normalization implies that when E does not occur, 
there may still be nontrivial evolution. In the continu- 
ous time paradigm, where one Kraus operator is given 



by Ey/dt, the transform (|2.9|) tells us that the no-jump 
normalization Kraus operator is given by 

r^o = 1 - ^E^Edt - ^(e""^c - e"l'c^)dt - iHdt. (4.5) 

Now we choose the driving Hamiltonian 

H ^^DS+'-^ie~"f'c-e''f'c^). (4.6) 

This is a Hermitian operator because of H4.4I) . Then the 
total evolution due to Qq is just the identity, apart from 
a term proportional to D{1 — S), which annihilates the 
codespace. Thus for a state initially in the codespace, the 
condition (|4.4() suffices for correction of both the jump 
and no- jump evolution. 

A nice generalization may now be found for a set 
{cj} of errors such that Cj [with associated operator 
Dj as defined in (|4.3|) ] acts on the jth qubit alone. 
Since Dj is traceless, it is always possible to find some 
other Hermitian traceless one-qubit operator Sj such that 
{sj, Dj} = 0. Then we may pick the stabilizer group by 
choosing the single stabilizer generator 

= si (g) . . . (g) s„ (4.7) 

so that the stabilizer group is {1,5'}. As noted in 
Sec. Ill A1 this is not strictly a stabilizer group, as S may 
not be in the Pauli group, but this does not change the 
analysis. Choosing H according to this S such that 

^ = E ^D^S + - e*c]) (4.8) 

j 

will, by our analysis above, provide a total evolution that 
protects the codespace, and the errors will be correctable; 
furthermore, this codespace encodes n — 1 qubits in n. 

Note that we can now easily understand the n-qubit 
jump process error of spontaneous emission considered in 
SecHni Here, 7 = 0, 5 = X®", and Dj = 2k jZj. Thus 
(|4.4|) is satisfied, and the Hamiltonian (|3.13|l is derived 
directly from (|4.8|l . 

Moreover, one is not restricted to the case of one sta- 
bilizer; it is possible to choose a different Sj for each 
individual error Cj. For example, for the spontaneous 
emission errors Cj = Xj — iYj we could choose Sj as differ- 
ent stabilizers of the five-qubit code. This choice is easily 
made, as the usual generators of the five-qubit code are 
{XZZXI,IXZZX,XIXZZ,ZXIXZ} 2]. For each 
qubit j, we may pick a stabilizer Sj from this set which 
acts as X on that qubit, and X anticommutes with 
Dj = Zj. This procedure would be useful in a system 
where spontaneous emission is the dominant error pro- 
cess; it would have the virtue of both correcting spon- 
taneous emission errors by means of feedback as well as 
correcting other (rarer) errors by using canonical error 
correction in addition. 

We note that the work in this section can very eas- 
ily be modified to generalize the results of That 
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work has the same error model as ours: known jumps oc- 
curing on separate qubits so that the time and location 
of each jump is known; but [23| postulates fast unitary 
pulses instead of a driving Hamiltonian. Their scheme 
for spontaneous emission depends on applying the uni- 
tary X®" at intervals Tc/2 that are small compared to 
the rate of spontaneous emission jumps. They show that 
after a full Tc period, the no-jump evolution becomes 

JJ ^ ^-iTaHa/2j^®n^-iTaH^/2j^(g)n ^ g-^-/2 X)" i 1 

(4.9) 

Thus the application of these pulses acts, as does our 
driving Hamiltonian, to correct the no-jump evolution. 
The generalization from spontaneous emission to general 
jump operator cj for their case is simple; the code is the 
same as in the above one-stabilizer protocol, with single 
stabilizer equal to (|4.7|l . The fast unitary pulses are in 
this case also simply equal to (|4.7|l . 



B. Diffusive unraveling 

The case of white- noise feedback, where 7 ^ oo, is 
easily treated by recalling the master equation H2.17|l for 
white-noise measurement and feedback. It is clear that 
the first term in 12.17|l can be eliminated by choosing the 
constant driving Hamiltonian 



H = 



(4.10) 



which is automatically Hermitian. The problem then 
becomes choosing a feedback Hamiltonian F such that 
c — iF annihilates the codespace. The choice for F can 
be made simply by noting that if the codespace is stabi- 
lized by some stabilizer S, we can choose 



F ^ B~ iAS. 



(4.11) 



Now, note that the decoherence superoperator V acts 
such that 



V[xI + L]p = V[L]. 



(4.12) 



Then we know that V[c - iF] = V[xl + A{I - S)] anni- 
hilates the codespace. 

The only caveat is that F is a Hamiltonian and there- 
fore must be Hermitian. Then the choice (|4.11|) for F is 
only possible if the anticommutator of S and A is zero: 



{S, A} = 0. 



(4.13) 



Therefore, if we are given the measurement operator 
g-i0g = ^ -f ^ -I- iB, we must choose a code with some 
stabilizer such that condition (|4.13(l applies; then it is 
possible to find a feedback and a driving Hamiltonian 
such that the total evolution protects the codespace. 

At first glance, it may seem odd that the condition for 
feedback does not depend at all upon B. This indepen- 
dence has to do with the measurement unraveling: the 
diffusive measurement record (|2.12|l depends only upon 
2iA + x)- 



C. Diffusion as the limit of jumps 

It is instructive to show that the diffusive feedback 
process can be derived by taking the limit of a jump 
feedback process using the transformation (|2.9|l . This 
takes several steps, and we use the treatment in js^ as 
a guide. But to begin, note that the condition (|4.13|) 
follows by considering Eq. (|4.4|l in the limit 7 ^ 00, as 
the leading order term in D is proportional to A. 

Consider the jump unraveling picture with jump oper- 
ator C-I-7 for 7 large (but not infinite). Recall that in the 
error-correction picture given in Sec. IIIII we postulated 
a feedback Hamiltonian {dN/dt)V that produces a uni- 
tary correction e~*^ that acts instantaneously after the 
jump. In addition we will postulate a driving Hamilto- 
nian K that acts when no jump happens. In this picture, 
we will show that given the condition H4.13II . it is possi- 
ble to find asymptotic expressions for V and K so that 
the deterministic equation for the system preserves the 
stabilizer codespace. Finally, we will show that taking 
the limit 7 — > cx) leads to the expression for the feedback 
and driving Hamiltonians H4.10|l and H4.11() . 

Let us consider the measurement operators for the un- 
raveling with large 7 and H — 0. Following (|2.9|l these 
are 

^0 = 1- y[c7-c^7+(c + 7)^(c + 7)], (4.14) 

where we have assumed for simplicity that 7 is real. Now, 
including the feedback and driving Hamiltonians modifies 
these to 



n[{dt) 

%{dt) 



dte-'"^ {c + 



-iKdt 



Qoidt) 



1 - iKdt 



dt 



{&c + 2-fc + Y)- (4-15) 



Following Ref. [33|, expand V in terms of I/7 to sec- 
ond order: V — ^^1/7+^2/7^ where the Vi are Hermitian. 
Then expanding the exponential in (|4.15|l we get to sec- 
ond order 



= ^/di-1 



1 



7 



7 



- + -{A- 

7 7 



_ivl 

2 72 

iB - iVi 



{A + iB + 7) 



^iV^y2~iV2-t{A + iB)Vi) 



(4.16) 



A reasonable choice for Vi, by analogy to H4.11|l . is 
B - iAS. Following we al so use f^HT^ and 
to choose V2 and K; note that (|4.11|) is exactly the ex- 
pression we would expect for K from l|4.6|l in the limit as 
7 is taken to infinity. We will proceed to show that the 
choice fir V and K, 



Vi 



B-iAS 



(4.17) 
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V2 = ~{c^F + Fc)/2 (4.18) 
K = ~--f{B~iAS), (4.19) 

leads to the correct evolution to second order in 7. 
Now, the deterministic evolution is given by 



dp = fl'oP^ + ilipfi'i - p. 



(4.20) 



Substituting (|^?T3|) - (|1?TO|) into (|On|) to second order 
in 7, after some algebra, gives the deterministic jump 
equation 



dp = V[A{1-S)]p 



(4.21) 



which of course acts as zero on the codespace. 

Now we will show that taking the limit as 7 — )■ 00 leads 
to the feedback operators given in H4.10|l and H4.11|l . We 
saw in H4.17|l and (|4.18(l that the feedback Hamiltonian 
needed to imdo the effect of the jump operator c + 7 was 
just 



dN{t) (B-iAS c^F + Fc 



dt \ J 2 
Keeping terms of two orders in 7 gives 

ct F _|_ Fr 

Hf, = jiB-tAS)-'-^^ 
dN{t) - "f^dt 



(4.22) 



+ - 



jdt 



{B-iAS). 



(4.23) 



The last term just becomes the current Q{t) as 7 ap- 
proaches infinity, as in equation H2.11|l . Furthermore, we 
have not yet added in the driving Hamiltonian to the 
expression for the feedback. Doing so yields 

c''F + Fc 

= Q{t){B-iAS)- '' ^ ' (4.24) 

which is just what we obtained in the previous section. 
Thus we can see that this continuous current feedback 
can be thought of as an appropriate limit of a jump plus 
unitary correction process. 



D. Imperfect detection 

These results for feedback were obtained by assuming 
unit efficiency, i.e., perfect detection. Realistically, of 
course, the efficiency r\ will be less than unity. This re- 
sults in extra terms in the feedback master equations we 
have derived In the jump case, the extra term is 



p^(\-^)Y.{c,pc]-U,c,pc]u]) 



In the diffusion case it is 

1 — 77 
P = 



(4.25) 



(4.26) 



In both cases this results in exponential decay of coher- 
ence in the codespace. This is because the error cor- 
rection protocol here relies absolutely upon detecting 
the error when it occurs. If the error is missed (jump 
case), or imperfectly known (diffusion case), then it can- 
not be corrected. This behavior is, of course, a prop- 
erty of any continuous-time error correction protocol that 
depends on correcting each error instantaneously (e.g., 

[Bill 111). 

On the other hand, such behavior for Markovian feed- 
back is in contrast to the state-estimation procedure used 
in . The latter is much more robust under non-unit ef- 
ficiency; indeed, given non-unit efficiency it still works to 
protect an unknown quantum state without exponential 
loss "33^. This difference in performance occurs because 
state-estimation is a function of the entire measurement 
record, not just instantaneous measurement results, and 
thence does not propagate errors to the same extent that 
a Markovian feedback system does. Thus we can see 
that there is a certain tradeoff. Our Markovian feedback 
scheme relies upon calculational simplicity, but at the 
expense of robustness. The state-estimation procedure, 
conversely, is designed to be robust, but at the cost of 
computational complexity. 



V. UNIVERSAL QUANTUM GATES 

Given a protected code subspace, one interesting ques- 
tion, as in |2ll |. is to investigate what kinds of unitary 
gates are possible on such a subspace. For universal 
quantum computation on the subspace — the ability to 
build up arbitrary unitary gates on fc qubits — it suffices 
to be able to perform arbitrary one-qubit gates for all 
k encoded qubits and a two-qubit entangling gate such 
as controUed-NOT for all encoded qubits /i, v. Indeed, 
as is noted in 21], it is enough to be able to perform 
the Hamiltonians X^,Z^, and X^X^ for all p,v |34| . 
We will demonstrate that performing these Hamiltonians 
with our protocol is possible for the spontaneous emission 
scheme given in Sec. lIIII and then we will show how that 
construction generalizes for an arbitrary jump operator. 

Recall that the example in Sec. Illll has single stabilizer 
X®" and encodes n—1 logical qubits in n physical qubits. 
To find the 2(n — 1) encoded operations, we must find 
operators that together with the stabilizer generate the 
normalizer of X®" In addition, if these operators are 
to act as encoded X and Z operations, they must satisfy 
the usual commutation relations for these operators: 



{X^,Z^} = 

[Xf,,Z,y] = 0,Py^ I' 

[X^,X,] = [Z^,Z,]^0. 



(5.1) 



Operators satisfying these constraints are easily found 
for this code: 



^ _ j®^— 1 j^jigin — /i 
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X^X^ = (5 2) 

where we assume l<fi<i'<n — I. If we apply a 
Hamiltonian iJenc given by any linear combination of the 
operators in H5.2I) , the resulting evolution is encapsulated 
in the expression for fla, from H3.6(l : 

no = {l-^Kjdt)l 

3 

- ^ HjZj{l - X'^")dt - iHencdt. (5.3) 

j 

As the first term is proportional to the identity and the 
second term acts as zero on the codespace, the effective 
evolution is given solely by -ffenc as long as the state 
remains in the codespace under that evolution. But be- 
cause the encoded operations are elements of the normal- 
izer, as we saw in Sec. Ill Al applying H^nc does not take 
the state out of the codespace. Furthermore, our protocol 
assumes that spontaneous emission jumps are corrected 
immediately and perfectly, so jumps during the gate op- 
eration will also not take the state out of the codespace. 
Thus we can perform universal quantum computation 
without having to worry about competing effects from 
the driving Hamiltonian. 

The generalization to the scheme given in Sec. II V Al to 
encode n — I logical qubits in n physical qubits is easily 
done. First we note that for the stabilizer S given in the 
general scheme, we know that 

S = UX'^^Ul (5.4) 

for some unitary U = ^"^i Ui, so the encoded operations 
for that code are similarly given by 

x^, = i®^'-^Uf,xuli®''-t' 
x^x^ = /®^-ic/^x[/t/®'^-^-if7^xt/J;/®"-^(5.5) 

Now, from H4.6II we can see that the generalization of 
(lOl is 

1^0 = (1 - fdt)l -Y,9jDj{l - S)dt - iHencdt, (5.6) 
j 



for real numbers /, gi given by expanding the expression 
(|4.5|) . Again, since D{1 — S) annihilates the codespace, 
the effective evolution is given solely by Henc as long as 
the state remains in the codespace under that evolution. 
Again, Henc is made up of normalizer elements, which 
do not take the state out of the codespace; and again 
jumps that occur while the gate is being applied are 
immediately corrected and thus do not affect the gate 
operation. Therefore, universal quantum computation is 
possible under our general scheme. 

VI. CONCLUSION 

We have shown that it is possible to understand a par- 
ticular variant of quantum control as quantum error cor- 
rection. This method is very general in that it can correct 
any single qubit detected errors, while requiring only n 
physical qubits to encode n — I logical qubits. As a par- 
ticular example, we have shown how to correct for spon- 
taneous emission evolution using feedback and a driv- 
ing Hamiltonian, which allows less redundancy than has 
previously been obtained. We have additionally shown 
that universal quantum computation is possible under 
our method. 

We expect that this work will provide a starting point 
for practically implementable feedback schemes to pro- 
tect unknown states. The fact that only two qubits are 
required for a demonstration should be particularly ap- 
pealing. We also expect a more complete theoretical de- 
velopment. Fruitful avenues for further research include 
applying notions of fault-tolerance to this sort of quan- 
tum control. 
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